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^ ; ABSTRACT 

, General solutions of the maser polarization problem are presented for arbitrary ab- 

sorption coefficients. The results are used to calculate polarization for masers permeated 
by magnetic fields with arbitrary values of xb , the ratio of Zeeman splitting to Doppler 
linewidth, and for anisotropic (m-dependent) pumping. In the case of magnetic fields, 
one solution describes the polarization for overlapping Zeeman components, xb < 1- 
The xb —> limit of this solution reproduces the linear polarization derived in previous 
studies, which were always conducted at this unphysical limit. Terms of higher order in 
xb have a negligible effect on the magnitude of q. However, these terms produce some 
major new results: (1) The solution is realized only when the Zeeman splitting is suffi- 
ciently large that xb > \fSo] J s , where Sq is the source function and J s is the saturation 
| intensity (pumping schemes typically have Sq/J s ~ 10 _5 -10 -8 ). When this condition 

' is met, the linear polarization requires J/ J s ;> xb, where J is the angle-averaged inten- 

sity. This condition generally requires considerable amplification, but is met long before 
saturation {J / J s > 1). (2) The linear polarization is accompanied by circular polariza- 
tion, proportional to xb- Because xb is proportional to the transition wavelength, the 
circular polarization of SiO masers should decrease with rotation quantum number, as 
observed. In the absence of theory for xb < 1, previous estimates of magnetic fields 
from detected maser circular polarization had to rely on conjectures in this case and 
generally need to be revised downward. The fields in SiO masers are ~ 2-10 G and were 
overestimated by a factor of 8. The OH maser regions around supergiants have fields of 
~ 0.1-0.5 mG, which were overestimated by factors of 10-100. The fields were properly 
estimated for OH/IR masers (<^ 0.1 mG) and H2O masers in star- forming regions (~ 
15-50 mG). (3) Spurious solutions that required stability analysis for their removal in 
all previous studies are never reproduced here; in particular, there are no stationary 
physical solutions for propagation at sin 2 9 < | , where 9 is the angle from the direction 
of the magnetic field, so such radiation is unpolarized. These spurious solutions can be 
identified as the xb = limits of non-physical solutions and they never arise at finite 
values oi xb, however small. (4) Allowed values of 9 are limited by bounds that depend 
both on Zeeman splitting and frequency shift from line center. At xb <J 10~ 3 , the al- 
lowed phase space region encompasses essentially all frequencies and sin 2 9 > | . As the 
field strength increases, the allowed angular region shrinks at a frequency-dependent 
rate, leading to contraction of the allowed spectral region. This can result in narrow 
maser features with linewidths smaller than the Doppler width and substantial circular 
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polarization in sources with xb ^ 0.1. When xb > 0.7, all frequencies and directions 
are prohibited for the stationary solution and the radiation is unpolarized. 

Another solution describes the polarization when the Zeeman components separate. 
This occurs at line center when xb > 1 and at one Doppler width when xb > 2. The 
solution is identical to that previously identified in the xb — ► oo limit, and applies to 
OH masers around HII regions. A significant new result involves the substantial dif- 
ferences between the tt- and u-components for most propagation directions, differences 
that persist into the saturated domain. Overall, HII/OH regions should display a pre- 
ponderance of c-components. The absence of any 7r-components in W3(OH) finds a 
simple explanation as maser action in a magnetic field aligned within sin 2 9 < | to the 
line of sight. 

Subject headings: masers, magnetic fields, polarization, radiative transfer 



1. INTRODUCTION 

Maser radiation is frequently polarized. Polarization generation is usually attributed to the 
presence of a magnetic field, which plays a dual role in this process. First, the magnetic field intro- 
duces a quantization axis, enabling a meaningful distinction between magnetic quantum numbers. 
This occurs when the field is sufficiently strong that the gyro-rate exceeds all other relevant micro- 
scopic rates. Once this condition is met, which is virtually always the case in astronomical masers, 
further increases in field strength are irrelevant; as far as this aspect of the problem is concerned, 
the field strength can be ignored altogether. 

Second, the magnetic field shifts the energies of magnetic sub-levels by the Zeeman effect, 
splitting each line into components with different Am. The relevant dimensionless parameter for 
the significance of the Zeeman splitting Avb is its ratio to the Doppler linewidth Ai/p, 

In the last equality g is the Lande factor with respect to the Bohr magneton, A is the transition 
wavelength in cm, B is the field strength in Gauss and Av^ is the Doppler width in km s _1 . Because 
the Zeeman shifts are proportional to the field strength, the absorption coefficients for the various 
Zeeman components vary with xb, thus this aspect of the problem can be expected to introduce a 
direct dependence of the polarization on xb- 

A general theory for the dependence of maser polarization on xb has not been formulated thus 
far. Instead, ever since the seminal work by Goldreich, Keeley & Kwan (1973; GKK hereafter) all 
maser polarization studies, including the first two in this series (Elitzur 1991, 1993; papers I and II 
hereafter), were performed in one of two limits — either ^ ooor^B ^ 0. Absorption coefficients 
were derived from rate equations in which these limits were assumed beforehand, not by taking the 
appropriate limits of general expressions for arbitrary xb- This procedure seems adequate when 
the Zeeman components are fully separated because the absorption coefficients drop exponentially 
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in the line wings. Once xb greatly exceeds unity, each component can be treated as an independent 
line, the dependence on xb can be expected to be insignificant and xb — ► oo is a proper, physical 
limit. 

The case of overlapping Zeeman components, xb <C 1, is considerably more problematic. This 
case was always treated in the limit xb = 0, which is of course unphysical — in the absence of 
a magnetic field there is no quantization axis and no polarization. In essence, current theory is 
predicated on the implicit assumption that, although xb is finite, it is sufficiently small that xb = 
can be used, but only in the expressions for the absorption coefficients (which are then derived from 
rate equations that assume this limit at the outset). It should be noted that, although unphysical, 
the xb — * limit of a complete theory can still be legitimately explored; the actual polarization 
varies with xb m some fashion and may well contain x^-independent terms that provide finite 
values in that limit. But such a limit can be considered only within the context of a general theory 
in which the x b variation is shown explicitly. The deficiency of current theory is that this variation 
is completely unknown because the limit xb = is assumed at the outset in the calculation of 
the absorption coefficients, making it an uncontrolled approximation. It is impossible to assess 
how small xb must be for any results to be valid and whether any of them may actually reflect 
singularities that arise only at xb = 0. 

All the theoretical tools for development of a general theory for maser polarization are available 
by now. Litvak (1975) has formulated the general theory for transport of line radiation and paper 
II has clarified the procedures introduced in GKK for identifying maser stationary polarizations. 
These ingredients are combined in §2 to produce the general solutions for maser polarization for 
arbitrary absorption coefficients. The subsequent two sections utilize these results for masers in 
magnetic fields — §3 presents the solution when the Zeeman components separate, §4 when they 
overlap. §5 presents the general solution for polarization in sources where the pumping is in- 
dependent with an arbitrary degree of anisotropy. All the solutions are tabulated for handy reference 
and the results are summarized and discussed in §6. This final section also contains comparison 
with observations and can be read on its own, independent of the rest of the paper. 



2. POLARIZATION SOLUTIONS 

A full description of the electromagnetic radiation field involves the 4-vector of its Stokes 
parameters S = (I ,Q,U,V). The general transfer equation for S in the case of line radiation was 
derived by Litvak (1975) and is listed in appendix A. When the source terms can be neglected, 
which is the case for maser radiation, the radiative transfer equation is 

f = R-S, (2-1) 

where the matrix R can be read off equation |A3| . Neglect of the source terms, further discussed 
in §4.2.1, is justified in essentially all cases of interest. Typically, the brightness temperature of 
an astronomical maser exceeds its excitation temperature by at least eight orders of magnitude, so 
only the properties of the self-amplified terms are relevant. 
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2.1. Stationary Polarizations 

The polarization is characterized by II = (q, u, v ), the 3- vector of normalized Stokes parameters 
q = Q/I, u = U/I, v = V/I (see figure 1). Configurations that obey dYl/dl = have stationary 
polarization that does not vary with propagation and can be expected to describe observed maser 
radiation; all other polarizations will evolve until they lock into a stationary state. The condition 
of stationary polarization is 

q> = u' = v' = 0, (2-2) 

where a prime denotes derivative along the path. Because q = Q/I, q' = is equivalent to 
Q'/Q = I' /I, with a similar relation for U and V. Therefore, stationary polarizations obey 

r _Q' _u> _v> _ 

where A is some (unknown) factor, common to all four Stokes parameters. So stationary polariza- 
tions obey dS/dl = AS and comparison with the radiative transfer equation shows that they are 
the eigenvectors of the matrix R, with A the corresponding eigenvalues. All four Stokes parameters 
of a stationary polarization state obey an equation identical to that for the intensity of a scalar 
maser, the widely used model that treats the radiation electric field as a scalar, with A its (single) 
absorption coefficient. With minor changes, the extensive theory derived for scalar masers can be 
carried over to polarized masers. In particular, the Stokes parameters of stationary solutions grow 
as exp(J Xdl). 



(A- 




= K t q + K C V 


(A- 


K m )q 


= H 


(A- 


K m )u 


= 


(A- 


K m )v 


= K c . 



2.1.1. The Solutions 

Because physical solutions must have 1^0, the eigenvalue equations, obtained from equation 
A3, can be divided by / to produce a set of equations for the normalized Stokes parameters: 

(2-4a) 
(2-4b) 
(2-4c) 
(2-4d) 

Now multiply eq. (b) by ki, eq. (d) by n c and combine. Inserting K\q + k c v from eq. (a) leads to 

(A - K m f = K 2 l +K 2 c = K 2 p , (2-5) 

so the eigenvalues are A = n m ± n p . When k v / 0, the solution is 

\ = K m + K p , q= — , n = 0, v= — , (2-6) 

Kp Kp 

which will be called a Type 1 Solution. Because k p > 0, the solution with A = K m — k p has a lower 
growth rate and can be dismissed. Since q 2 + v 2 = 1, type 1 solutions describe fully polarized 



All the explicit type 1 polarizations derived below have hip — Km and this choice of sign gives a solution that 
does not grow at all. 
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radiation. Note that although we obtained explicit expressions for the fractional polarizations, m 
and k c themselves can depend on q and v. In that case the polarization may not be fully specified, 
with the solution only providing a self-consistency relation for q and v. 

When K p = (i.e., both k\ and k c vanish), A = K m . All the eigenvalue equations are obeyed 
trivially and the polarization is not constrained by the eigenvalue problem. However, only certain 
polarization configurations can maintain K[ = k c = 0, if at all, so the polarization is determined in 
this case from the conditions imposed by the vanishing of the off-diagonal elements of the matrix 
R. That is, the solution is 

A = K m , m = k c = 0, (2-7) 

which will be called a Type Solution. Just as in the general case of type 1 polarization, the Stokes 
parameters are determined implicitly. Because the off-diagonal ki and k c determine all the possible 
differences among the three basic absorption coefficients K Am , their vanishing implies 

K + = KT = K° = K m = A. (2-8) 

Type polarization is obtained when the absorption coefficients of all three Am transitions are 
equal to each other. In the case of J = 1 — > transitions, discussed in appendix A, this implies 
equal populations for all the magnetic sub-states of the upper level, a point first noted in paper I. 
Obviously, such a situation is impossible when only one Am transition is involved, as is the case for 
fully developed Zeeman patterns. Type solutions, which can be identified using either of the last 
two equations, are possible only when the Zeeman components overlap. Note that such a solution 
is not the k p — ► limit of a type 1 polarization; although the eigenvalue is properly reproduced, q 
and v are indeterminate in that limit. 



2.2. Miscellaneous Properties of the Solutions 

The polarization solutions can be further understood with the aid of 

k p = (ki,0,k c ), (2-9) 
a vector in (q, u, v )-space whose magnitude IS I f^p I — ^-p 

(see figure 1; this vector was denoted (3 in 

Litvak's paper); note that K p • II = K\q + k c v. With this vector, the radiative transfer equations 
can be written as 

^ = ( K m + K p • n)/ 

^ = Kp - (k p ■ n)n = [Uxk p }xu + (i -n 2 )K p . (2-10) 

The polarization is stationary (II' = 0) when K p vanishes, corresponding to a type solution. When 
Kp 7^ 0, the polarization is stationary only when II is along K p and has a magnitude of unity (full 
polarization), corresponding to a type 1 solution. This result has a simple interpretation in terms of 
the individual polarization modes II n , which make up the overall polarization through the ensemble 
average II = (II n ) (see paper II). Because |II n | 2 = 1, the radiative transfer of individual modes 
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involves pure rotation with the angular velocity Yl n XK p , which is different for different modes. The 
rotation comes to a halt only when II n overlaps k p . The overall polarization II settles into the 
stationary solution when all individual modes have rotated into overlap with n p , producing a fully 
polarized radiation field whose polarization is determined by evident from fig. 1. 

With this form, the radiative transfer equations admit the formal solution 

I = I e Tm+T *, n 2 = 1 - (1 - ng)e~ 2T % (2-11) 

where 

i r-i 

K m dl, tv = / II ■ Kpdl (2-12) 
la J la 

and where Iq and IIo are the intensity and polarization at some fiducial point Iq along the path. The 
vanishing of n p produces r n = and stationary polarization whose intensity increases exponentially 
with r m , i.e., a type solution (eq. |2-7|) . In addition, full polarization always maintains a fixed 
magnitude. It maintains also a fixed direction, becoming stationary, when II overlaps n p , leading 
to 7v = J tv p dl and the type 1 solution of equation p-6| . 



2.2.1. Scalar vs Polarized Masers 

The general formal solution of the radiative transfer problem helps illustrate a fundamental 
difference, first noted in Elitzur (1995), between the radiation fields of scalar and polarized masers. 
In the scalar maser model, different radiative modes correspond to different initial intensities Iq 
and the growth of each mode follows the formal solution 

I = h e T , (2-13) 

where r = J ndl and k is the (only) absorption coefficient, the equivalent of K m . In the unsaturated 
regime k is a constant independent of intensity, and the growth is exponential with pathlength. 
Therefore, the scalar unsaturated maser is a purely linear amplifier; its growth rate is independent 
of intensity, i.e., it is mode invariant. All the modes grow at the same rate and the mode distribution 
is unaffected by the amplification. This is the reason why the need to consider ensemble averages 
usually does not even arise in the scalar maser model. 

The inclusion of polarization changes the situation fundamentally. The formal solution for / 



is now given in eq. 2-11 and similar expressions can be listed for all other Stokes parameters. As 
with the scalar maser, when radiative interactions are negligible in the unsaturated domain the 
absorption coefficients are constant and their integrals increase linearly with pathlength, resulting 
in exponential amplification. However, even though the absorption coefficients are independent of 
the Stokes parameters in this case, the optical depth r n still does depend on them. Because of 
the explicit dependence of T n on II, the growth of the Stokes parameters depends on the Stokes 
parameters themselves, in spite of the fact that the population inversions are unaffected by radiative 
interactions. In order to get a quantity whose growth is mode invariant, similar to the intensity of 
a scalar maser, the off-diagonal optical depth r n must be eliminated. This is easily accomplished 
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by combining the expressions for 7 and II 2 in eq. [2-ll| to produce 



7(1 - n 2 ) 1/2 = 7 (1 - III) 1 ' 2 e Tm , (2-14) 



the analog of eq. 2-12 for scalar masers. In the unsaturated domain, r m again is the same for all 
modes. However, now the invariant growth does not apply to any single Stokes parameter, instead 
it applies to 7 2 (1 — II 2 ) = I 2 — Q 2 — U 2 — V 2 , a quadratic form involving all four Stokes parameters.! 

Since the growth rates of the Stokes parameters depend on the Stokes parameters themselves, 
the unsaturated polarized maser is a non-linear amplifier. As a result, a calculation of the evolution 
of ensemble-averaged Stokes parameters, rather than merely deriving their stationary values as done 
here, invariably requires numerical simulations that properly evolve an ensemble of waves. Such 
a simulation, frequently performed in plasma studies, will not be attempted here. It would be 
similar to a demonstration of the approach to Maxwellian of a particle velocity distribution or the 
approach to Planckian of a photon distribution. While such demonstrations are valuable, they 
are considerably more difficult than the mere derivation of either equilibrium distribution, the 
equivalent of the stationary solution of the maser polarization problem derived here. 

Stationary polarization requires that the radiation maintain the same normalized Stokes pa- 
rameters during propagation through the source, even as the Stokes parameters themselves grow by 
maser amplification. However, the solutions for such polarizations are determined by the vector n p , 
and in general K p can be expected to vary both in magnitude and direction, reflecting variations of 
the level populations as a result of interaction with the growing maser radiation. A solution is not 
guaranteed; indeed, stationary polarization may be altogether impossible in any particular region 
of phase space, in which case the radiation is unpolarized. We proceed now to derive the stationary 
polarizations, wherever they exist, for various magnetic field strengths. 



3. FULLY RESOLVED ZEEMAN PATTERN 

When the Zeeman pattern is fully developed, only one absorption coefficient is different from 
zero at the frequency of each component. Thus the condition cannot be fulfilled 

and type polarization is impossible; only type 1 polarization can grow in this case. In addition, 
the direction of the vector K p (though not its magnitude) is constant, unaffected by the variation 
of the level populations with saturation. 

At the frequency of a 7r-component, only kP does not vanish. Therefore, k c — while K m — 
K p = —K[ = ^K°sin 2 6. From eq. |2-6| , the polarization solution is 

K n = k° sin 2 9, q = -l, v = 0, (3-1) 



2 This result can also be obtained directly from the radiative transfer equation 



A3 which leads to 



-(I 2 - Q 2 - U 2 - V 2 ) = 2K m (I 2 -Q 2 -U 2 - V 2 ). 
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where k 71 (= K m + k p ) is the eigenvalue, i.e., the maser growth rate. In the case of cr-components 
only k + (k~) does not vanish, so Ki = l^sin 2 ^, k c = ihn^cosO and K m = k v = \k^(1 + cos 2 6). 
From eq. |2-6| , the polarization solution is now 

1+/ 9 „s sin 2 6* 2 cos # 

K ff = l« ± (l + COS 2 e), </= — ^ = ±T- 2fl, 3-2 

z 1 + cos 1 + cos 6 

where kF is the appropriate growth rate. These polarizations are identical to those derived in GKK 
for fully resolved Zeeman patterns. Since the detailed expressions for the absorption coefficients K Am 
were not even specified, these solutions apply to all spins and degrees of saturation, in agreement 
with the conclusions of paper I. 

Explicit expressions for k 71 and K a can be derived only in the context of a specific maser 
model. Appendix A provides the general solution for the absorption coefficients K^ m of a spin 



1 — > transition (equation A7). When the Zeeman components are fully separated, ground state 
particles with a given velocity interact only with radiation shifted by the same amount from each 
of the three line centers; for example, particles with zero velocity interact only with the line center 



radiation of each Zeeman component. Therefore, each J/\ m and Kq 171 in eq. A7 must be considered 



at the same frequency shift from the line center of the (l,m) — ► transition. If pumping of the 
upper levels is m-independent, i.e., P\ m = -Pi, then Kq 171 = kq oc (Pi — Po)/T and the absorption 
coefficients are 

KAm = 1 I f K °r TT E ' where E = f 1 + E ■ ( 3 - 3 ) 

1 + JAmJAm/ J.s \ Am J 

The common factor E varies from 1 when all three components are unsaturated {e& m = 0) to 
j when all are saturated (eAm = !)• Inserting the polarization from the appropriate solution 



(equations 3^1 and p-2|) , the absorption coefficients of the Zeeman components become 



K sin 2 6> kqxUi + cos 2 9) 

« = % 7— E > K<T = TT^ o \ , E, (3-4) 

1 + sin 2 9J n /J s 1 + |(1 + cos 2 9)J a /Js 

where J n and J a are the corresponding angle-averaged intensities. These results provide reasonable 
approximations for the absorption coefficients at frequency x when xb > 1+x because the matching 
spectral segments of neighboring lines are then separated by at least a full Doppler width. Therefore, 
these results apply at each line center when xb > 1, and as xb increases their applicability region 
spreads to additional frequencies, further removed from line center. When xb ^ 2, these absorption 
coefficients should describe adequately the entire portion of interest of each line. 

Apart from the factor E, for propagation along the axis the absorption coefficients of the cr- 
components are identical to those of a scalar maser with the same pumping scheme. The same 
applies to the 7r-component, but for propagation perpendicular to the axis. Moving to other direc- 
tions, in each case the unsaturated absorption coefficient is reduced by an angle-dependent factor. 
When both components are unsaturated, their intensities grow according to 

In = SqX exp(nol sin 2 9), l a = Sq exp(^K /) x ex.p(^Kol cos 2 9), (3-5) 
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where Sq is the source function. These angular distributions are highly anisotropic; with kqI = 10, 
corresponding to unsaturated maser operation in typical pumping schemes, I a drops to 1/e of its 
peak at 9 = 27° and the same holds for I T at 90° — 9 = 18°. The anisotropy introduces a high degree 
of asymmetry between the it- and cr-components, which have equal intensities for the same growth 
length at sin 2 9 = |, i.e., 9 = 55°. At < 9 < 55° the cr-components have a higher intensity, at 
55° < 9 < 90° the 7r-component is stronger. The regions where each Zeeman component dominates 
have comparable sizes, but the enhancement of the cr-components in their region of dominance 
is more prominent. Moving away from 9 = 55°, the ratio n 71 / 1 K a reaches a maximum of 2 for 
propagation at 9 = 90° while the inverse ratio n a / becomes 2 already at 9 = 39°, increasing 
further without bound as the magnetic axis is approached. Therefore, an overall preponderance 
of cr-components can be expected among all unsaturated maser sources with comparable physical 
conditions. 

When both components saturate, the ^-dependence of their absorption coefficients disappears 
and they both obey kJ = \kqJ s oc Pi — Pq. Saturated regions have equal production rates of a- 
and 7r-photons independent of inclination to the magnetic axis. However, in most cases the lengths 
of such regions will be ^-dependent and different for the two components. For each transition, 
the unsaturated absorption coefficient is reduced by a 6>-dependent factor and the corresponding 
saturation intensity is increased by the same factor. Both effects cause the optical depth that brings 
saturation to increase, lengthening the unsaturated region. If t s denotes the optical depth required 
for the saturation of a scalar maser in the linear geometry, the corresponding saturation optical 
depths of the Zeeman components are 

_ n T s -lnsin 2 (j r s - In ±(1 + cos 2 0) 



sin 2 9 s ±(l + cos 2 0) 



These results follow from equation |3-4| and general expressions for linear masers (Elitzur 1990). 
They show that saturation, and ^-independent maser production, generally requires pathlength 
longer than t s . But whereas the saturation optical depth of cr-components can only be as large as 
~ 2t s , for the 7r-component it can become much larger, diverging when 9 — > (reflecting the fact 
that this component does not grow for propagation along the axis). Once the maser length in a 
given direction exceeds that required for saturation of the weak component, the intensities of both 
components continue to grow at the common rate jKoJ s per unit length, maintaining a constant 
difference 

J a -J n = ij s (rJ-r s CT ). (3-7) 

The indirect coupling through the common factor E (eq. [3-3| ) introduces an additional element 
of asymmetry into the maser growth pattern. In a region where one Zeeman component is saturated 
but the other is not, E is smaller than its unsaturated limit of 1 but larger than its saturated limit 
of j. Therefore, the growth rate of the stronger, saturated component is larger than what it would 
be were all components saturated, at the expense of the weaker component whose unsaturated 
growth rate is suppressed. When only the cr-components are saturated, Kq is reduced by a factor 
of 3, and when only the 7r-component is saturated, Kq is reduced by a factor of 2. This suppression 
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of the weak Zeeman component by the strong one should add to the dominance of cr-components 
in maser sources with xb > 1- 



4. OVERLAPPING ZEEMAN COMPONENTS 

When xb < 1, the overlap of Zeeman components engulfs most of the line cores. In the maser 
model discussed in appendix A, ground state particles at a given velocity can interact with different 



Zeeman components at the same frequency so all quantities in equations A7 for the absorption 
coefficients must now be considered at a given frequency; in particular, J& m = J, the angle-averaged 
intensity at that frequency. Because of the Zeeman splitting, a given frequency corresponds to 
different shifts from the line centers of the three Zeeman components. The (1, m) — * transition is 
centered on the frequency v m = vq — ttiAvb and the argument of the corresponding Doppler profile 
is (u — v m )/AvD = x + mxB, where x = [y — vq)/Aud is the dimensionless frequency shift from 
the center of the unperturbed line. Assuming m-independent pumping, the unsaturated absorption 
coefficients of the three transitions are equal at their respective line-centers, KQ Lm (is m ) = kq oc 
(Pi - P )/r, so 

K fr m (x) = k exp[-(x + x B Am) 2 }. (4-1) 

In spite of the m-independent pumping, the ratios of absorption coefficients vary across the line 
because of the Zeeman shifts. These variations are conveniently expressed in terms of the ratios 

Ri = K o + *o = e -x% cos h(2xx B ) ~ l-x 2 B (l-2x 2 ) 

R c = K ° = - e - x B S mh(2xx B ) ~ -2xx B - (4-2) 

2k 

In each case, the second equality holds for any xb while the last relation provides the appropriate 
small-x^ limit, a result valid to second order in x# and all frequency shifts with 2xx# < 1. The 
constraints imposed on k^™ by the inequalities of eq. A12| produce 1 1 — Pi | < \ , or 



x| < 



2 1 - 2x 



(4-3) 



At any frequency x, the Zeeman splitting must obey this bound to ensure that all three transitions 
remain inverted when they are saturated. In particular, fulfillment of this constraint over the entire 
central portion of the line, |x| < 1, restricts the Zeeman splitting to xb < l/v2- 

When the three unsaturated absorption coefficients are expressed in terms of k{] and the ratios 
Pi and P c , the explicit results for the absorption coefficients (eq. |A9| ) become K^ m = KqK^ 171 /D, 
where K^ m are dimensionless quadratic forms in J/J s - With these expressions, K m , ki and k c 
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take similar forms whose numerators are 

K m = i[sin 2 6» + i?i(l + cos 2 6») 

+ (J/J S ){(1 + cos 2 0)[/i(2 - 3i?i) + 2R 1 - 1] + 2sin 2 0/i(2 - - wR c (5 + cos 2 0)} 
+ (J/J S ) 2 {(1 + cos 2 9)f [f 1 (2R 1 - 1) - 4wR c ] + sin 2 0(2i?x - 3)(w 2 - f 2 sin 2 9)}} 
K x = ±sin 2 #[i?i-l 

+ (J/J s )[2i?!-1- + 2)- wi? c ] 

+ (J/J s ) 2 {/o[/i(2i?i - 1) - AwR c ] - (f 2 - w 2 )(3 - 2R,)}} 
K c = cos 9[R C + (J /J S )[R C (2- f 1 )- W R 1 ] + (J/J s ) 2 f [4f 1 R c -w(2R 1 -l))}. 

(4-4) 

In each of these relations the J-independent term provides the J = (unsaturated) limit of the 
corresponding absorption coefficient (since the J-independent term in D is 1). There is a funda- 
mental difference between these terms for the diagonal and off-diagonal absorption coefficients when 
xb < 1- While K m (J=0) ~ 1 in this case, |-K";(0)l an d |-PsT c (0) | are proportional to \R\ — 1\ ~ x\ 
and \R C \ ~ xb, respectively, and thus are much smaller than unity. As a result, although the 
J-independent term provides an adequate approximation for n m in the entire unsaturated domain 
(J/J s < 1), ki and k c are dominated by J-dependent terms long before saturation (J/J s = 1); that 
is, the dependence on J of the off-diagonal absorption coefficients becomes significant at intensities 
well inside the unsaturated domain, either J/ J s ~ x\ or J/J s ~ ig. The reason for this difference 
is simple. While the J-dependence of K m reflects only the difference between the mean populations 
of the upper and lower levels, the corresponding dependence of k\ and k c reflects also population 
differences among the magnetic sub-states of the upper level. And those are affected by radiative 
interactions long before saturation reduces the overall population inversion. 



4.1. Type 1 Polarization 

Type 1 solutions can be derived from the equation qK c (q,v, J/J s ) = vKi(q,v,J/J s ), where 
v is replaced everywhere by \/l — q 2 (see eq. 2-6). Since different powers of J/J s have different 
coefficients in K\ and K c , the vector K p rotates in the q—v plane during the growth to saturation. 
Different results are obtained in the unsaturated (J/J s = 0) and saturated (J/J s — ► oo) limits, and 
a single type 1 solution for all intensities does not exist; if a proper physical solution is to exist, it 
must be able to evolve from the former to the latter. 



4- 1-1. The J = Limit 



In this limit, the absorption coefficients of eq. 4-4 are applicable for every value of xb- The 
reason is that, in the absence of radiative interactions, the rate equations have an identical functional 
form for all Zeeman splittings. In this limit K\ and K c are independent of q and v and the 
polarization solution can be obtained by inserting the appropriate terms from eq. 4-4 directly into 
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eq. 2-6, yielding 

(gi - l)sin 2 6> 2R C cos g 



- l) 2 sin 4 + 4#2 cos 2 _ !)2 sin 4 + AR 2 cos 2 



(4-5) 



This result provides the type 1 polarization at J = for any Zeeman splitting and can be used to 
follow the variation of polarization with magnetic field strength at any frequency and propagation 
direction. In particular, when xb > 1 these expressions should reproduce the solutions for the fully 
resolved Zeeman case, which are J- independent, at the three frequencies v m . This is indeed the case. 



Comparison with eqs. |3-1| and |3-2j shows that the 7r-polarization is produced when R\ = R c = 0, the 
<T-polarizations when both R\ and R c greatly exceed unity with their ratio obeying R\/R c = ±1. 
Indeed, these are the respective behaviors of R\ and R c when xb > 1- Consider first the ir- 
component. In this case, both R\ and R c vanish across the entire line (|x| <J 1) for large xb, as can 
be seen from eq. 4-2: at x = 0, R c vanishes for any xb and R± = exp(— x B ), vanishing rapidly when 



xb > 1; and at \x\ = 1, both i?i and R c vanish faster than exp(— xb) when xb > 2. Therefore, 



the 7r-polarization is properly reproduced by eq. 4-5 and is established at the line wings at slightly 



larger Zeeman splitting than at their centers, in agreement with the conclusion of the previous 
section. And at the frequencies of the a-components, 

x = tx B : Ri = \e x l (l + e~<) , R c = ±\e x l (l - e~<) . (4-6) 

Therefore, both ratios increase rapidly and obey R\/R c = ±1 when xb > 1, properly reproducing 



the o"-polarizations of eq. |3-2 



Proceed now to < 1- The polarization properties of the solution are quite different in this 



regime because both R\ — 1 and R c are small quantities in this case. With the aid of eq. 4-2 , 
expansion to second order in xb yields 

„ _ -zb(1 ~ 2x 2 )sin 2 g >( -4x cos 9 ( 4 _ 7 ) 



x 2 B (l - 2x 2 ) 2 sin 4 6 + 16x 2 cos 2 6 \f - 2x 2 ) 2 sin 4 9 + 16x 2 cos 2 9 

When 4|x| cos 9 > xb sin 2 9, the polarization is virtually purely circular, v ~ ±1, with the opposite 
sense in the two halves of the findl. The transition from one sense of circular polarization to 
the other occurs over the central frequency region A\x\ cos 9 < xb sin 2 9, where the polarization is 
purely linear, q ~ — 1. This fully polarized configuration requires precise phase relations among the 
three Zeeman components. Since the three transitions are pumped independently, this solution is 
unstable against perturbations, a conclusion confirmed by formal stability analysis. 

4-1-2. The Saturated Limit 

The polarization relevant for observations involves the saturated limit, J/J s S> 1, and can 
be obtained by retaining only the terms quadratic in J/J s of the equation qK c = vK[. From the 



3 By contrast, pure cr-components generate ^-dependent elliptical polarization (eq. p-3). 



-13- 



explicit form of the resulting 4-th order equation, an overall term q — qi factors out where 

sin 2 8 2 cos 8 , , , 

1 + cos z 1 + cos z 

This can be recognized as the <7-polarization in the fully resolved Zeeman case but with the opposite 
sign for q (see eq. 3-2), providing finite circular polarization that is independent of both xb and x. 



From symmetry this is impossible, so this is obviously a spurious solution. Indeed, there is no path 
to this polarization from the J = limit. 

After factoring out the term q — qj, the remainder is a cubic equation for q. All three roots of 
this equation involve the square root of a negative term so there are no physical solutions; type 1 
polarization is impossible now. In each case, the coefficient of the imaginary part is proportional to 
R c , which vanishes at xb = 0. Therefore, at this unphysical limit the imaginary parts are removed, 
producing the three type 1 solutions 

TT ■ „ = l + 3cos 2 fl „ = , 2(2 + 3 cos 2 8) 1 ' 2 

q 3(1 + 008*0)' 3(1 + cos 2 8) 

III : q = +1, v = 

IV : q = -1, v = 0. (4-9) 

Indeed, these solutions were identified in previous studies, which were always conducted at the 
xb = limit, and were shown to be unstable against perturbations. Solution II was found by GKK 
(their eq. 61), who showed it to be unstable for all 8. Solution III was found in paper II (eq. 2.3) 
and likewise was shown to be unstable for all 8. Solution IV was first derived by GKK (eq. 60), 
who found it to be unstable for sin 2 8 > ^ but stable for sin 2 8 < |. However, paper II showed this 
solution to be unstable also for sin 2 8 < | in the unsaturated domain, prohibiting its growth. The 
present general analysis, which is not restricted to any particular value of xb, shows that all three 
are simply the real parts of spurious solutions whose imaginary parts vanish only in the unphysical 
limit xb = 0. These polarizations are not realized at any finite value of xb, however small. 

In conclusion, there is no stationary type 1 solution that grows into the saturated regime when 
xb < 1. Just as type polarization is impossible when the Zeeman pattern is fully resolved, type 
1 polarization is impossible for overlapping Zeeman components. 



4.2. Type Polarization 

The type solution can be derived from K\ = K c = or, equivalently, from 
Either condition yields 

„. _ 8R r 



1 



(2Ri + l)sin 2 



3 _ 2 R X + 2(1 ~ R ' 



J/Js 



{2R X + l)cos0 V ^ 4J/J. 



1 + TTTT ■ (4-10) 



The corresponding eigenvalue, the common growth rate of the four Stokes parameters, is 

(2gi + lV3 „ jl- §4(1-23 



- ^ tTfJU, ■ (4 - n) 



-14- 



a ^-independent growth rate similar to that of a scalar maser with the same pumping scheme; note 
that K°(2iii + l)/3 = §04 + k + «g). 

Since the three absorption coefficients K Am are different from each other when J —* 0, obviously 
the type condition cannot be obeyed in that limit. A fully stationary polarization is impossible 
during early stages of maser growth, and this is the reason why both q and v of the formal solution 



contain J-dependent terms. However, as seen from eq. |4-4| , relative differences among absorption 
coefficients at J = are only of order at most and can be overcome when radiative interactions 
become important. Indeed, the J-dependent terms in eq. 4-10| vanish at high intensities and can be 



neglected in the expression for q when J/J s > 2|1 — Rx\ ~ 2x 2 B \l — 2x 2 \, for v when J/J s > |. Once 
these conditions are obeyed, the J-independent parts of the solution properly describe stationary 
polarization. Expanding the ratios R\ and R c to second order in xb, the stationary polarization 
solution is 



q = l ^r- 1 + f x 2 (1 - 2x 2 ) , v = - ^4 . (4-12) 

3 sin 2 9 L 35 J 3cos# v ' 

The inverse dependence of v on cos 9 reflects the fact that the circular polarization enters into the 

problem only as w = \vcos9 (see eq. |A5|) . 

Physical polarizations must always obey q 2 + v 2 < 1. Applied to the stationary solution, this 

inequality translates into a quadratic equation for sin 2 6, placing bounds on the allowed propagation 

directions 

|[1 + §x|(l + 2x 2 )] < sin 2 < 1 - 32x 2 B x 2 . (4-13) 

Figure 2 displays these bounds by plotting the largest and smallest allowed angles as functions of 
xb for various frequency shifts. At each of the marked values of x, physical solutions are possible 
only inside the 6 — xb region bounded by the corresponding curves. At xb = these bounds yield 
only the constraint sin 2 9 > I, i.e., 9 > 35°, familiar from the GKK study. However, at any physical 
value of xb the upper bound is also meaningful, and propagation perpendicular to the magnetic 
axis is forbidden (except for line center). As xb increases, the lower bound on sin 2 9 increases 
and the upper bound decreases, causing the allowed region of propagation directions to shrink. 
Finally, a physical solution is no longer possible when the two bounds coincide, and the radiation 
is unpolarized across the entire line for all propagation directions. 

As discussed before, strictly stationary polarization is impossible during the early stages of 



maser growth. From the formal solution of the radiative transfer problem (eq. |2-11 ), the amount 



by which any polarization varies during maser growth is controlled by the magnitude of the off- 
diagonal optical depth 7V for that polarization. A polarization II is approximately stationary if 



for that polarization is small. For the polarization solution of eq. 4-12 , T n is related to the diagonal 
optical depth r m through 

Tn = \T m x 2 B [l + 30x 2 + 3(x 2 - \) sin 2 9]. (4-14) 

This result is valid to order x B in the entire unsaturated domain, J/J s < 1, the region where 
optical depths grow linearly with pathlength. It shows that for radiation polarized according to the 
polarization solution, r n is much smaller than r m , the optical depth that controls the exponential 



-15- 



growth of the Stokes parameters. To a good degree of approximation, the polarization can be 
considered stationary if \r n \ remains less than unity for the entire unsaturated growth phase. This 
phase is characterized by r m = r s , where r s is the optical depth that brings saturation, thus the 
polarization solution is valid so long as \r n \ < 1 for r m = r s , or 



x 2 B \l + 30x 2 + 3(x 2 - \) sin 2 0\ £ |-. (4-15) 



This constraint ensures that the solution is approximately stationary also for the early growth 
phase. Figure 2 displays the bounds imposed by this constraint, too, for t s = 15, which is typical of 
astronomical masers' pumping schemes. As is evident from this figure, which plots also the bound 



of eq. 4-3, the primary bound at all frequencies is provided by eq. 4-13 , with the other two bounds 
providing further minor restrictions at various frequencies. If we require propagation up to 9 ~ 80° 
at all frequencies \x\ < 1, then the Zeeman splitting allowed for type polarization is limited to 
x B & 0.03. 



4-2.1. Source Terms and Onset of the Solution 

The polarization solution cannot be established before the J-dependent terms start dominating 
the behavior of K\. This occurs when the term linear in J becomes equal to the J- independent 



term, i.e., J/J s — x B (see eq. 4-4). However, the term that eventually controls the behavior of K\ 
once the maser saturates is the one proportional to (J/J s ) 2 . Indeed, the stationary polarization 
solution reflects the coefficient of this term and thus cannot be established before it too becomes 
comparable to the J-independent term, i.e., (J/J s ) 2 ~ x 2 B . Thus the solution requires J/J s ^> xb- 
Source terms were ignored thus far. Pumping of astronomical masers typically has So/Js ~ 
10 -5 -10~ 8 , thus the maser intensity greatly exceeds the source function So at early stages of 
growth, long before saturation. However, the polarization solution reflects the behavior of the off- 
diagonal terms in the transfer of self-amplified radiation while the source terms involve the diagonal 
absorption coefficient K m . Thus the onset of the linear polarization solution requires K[I > n m So, 
or, since K m ~ 1 in the entire unsaturated domain, k\I j ' J s > So/J s . In linear masers this condition 
is equivalent to KiJ/J s > Sq/J s , in three dimensions there is an additional factor of £1/Att, where Q 
is the beam solid angle. This condition must be obeyed when the J-dependent terms dominate the 
behavior of ki , at which time m is of order J/ J s . Thus the dominance of self-amplified radiation 
requires (J/J s ) 2 > So/J s , a condition first noted in paper II. For this condition to be obeyed when 
the quadratic term surpasses the J-independent term (J/J s ~ xb), the Zeeman splitting must obey 
x 2 B > So/Js! weaker fields will produce only unpolarized radiation. Combining these results yields 



J/Js ^x B > J So/Js ■ (4-16) 



For a given pumping scheme, these are the general bounds that must be obeyed by the Zeeman 
splitting and intensity to enable the polarization solution. 
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5. ANISOTROPIC PUMPING 

An alternative to magnetic fields in generating maser polarization is anisotropic interaction 
rates, where a quantization axis is introduced by the pumping process itself. Among the concrete 
examples proposed are collisions with an electron stream (Johnston 1967). Because the pumping is 
then m-dependent it directly introduces distinctions among magnetic sublevels even though they 
remain strictly degenerate. As long as symmetry is maintained around the and 
k c = R c = v = 0. Circular polarization is impossible under these circumstances. 

Consider the 1 — > maser model. The absorption coefficients are derived in the same manner 
as for overlapping Zeeman components but the ratio R\ is now frequency independent, given by 

Rl = Twv (5 - 1} 



Here P\\ denotes the common pump rate of the m = ±1 sublevels. From the bounds of eq. A12 



R\ must obey \ < R\ < | whenever all three transitions are radiatively coupled, to ensure that 
they remain inverted in the saturated domain. In addition, R\ = 1 is unphysical because an 
axis cannot be defined in that case. However, similar to xb = for magnetic fields, this value 
can be approached as a limit since arbitrarily small |1 — R\\ correspond to physical situations. 
The absorption coefficients again can be written as K Am = KqK m /D, and the corresponding 
expressions for D and K^ m are simpler than in the case of Zeeman splitting. The relevant analogs 



of eq. 4-4 are 

D = l + (2- fi)J/J s + 4fofi(J/Js) 2 

Ki = \ sin 2 6{R X - 1 + (J/ J s )[2Ri{l - /i) - 1 - /i]). (5-2) 



5.1. Type 1 Polarization 

Because k c = 0, n p = \m\ and from eq. |2-6| , the type 1 polarization is 

\ K i\ \ K i\ 

Self-consistency requires K\ to have the same sign as q. Since the J-independent term of K\ is 
proportional to R\ — 1 (eq. |5-2| ), unsaturated growth is possible for q = +1 polarization only when 
Ri > 1, for q = —1 when Ri < 1; the discontinuous transition between the two solutions occurs at 
the singular point R\ = 1. Furthermore, K\ must maintain its sign after the maser saturates. 

Consider first the q = +1 solution. In this case fi = \ and K\ remains positive after saturation 
only if Ri > I. Although the right-hand inequality of eq. A12 is violated in that case, the solution 
is still valid. That particular inequality reflects the requirement that the Am = transition remain 
inverted in the saturated domain. But this requirement does not apply here because q = +1 implies 
fo = 0, and the Am = transition does not couple radiatively. This polarization is generated purely 
by the |Am| = 1 transitions, indeed its growth rate is 

* = ^ + « = TTWj.> (5 - 4) 
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where Kq (= KqR\) is the unsaturated limit of re 1 . 

The q = — 1 polarization grows when h < R\ < 1 and produces /o = sin 2 9, fi = ^cos 2 #. 
The requirement that K\ remain negative after saturation does not further constrain R\, instead it 
restricts the propagation directions to 

sin 2 0< ^L^i. (5.5) 

This range in displayed in figure 3. At the singular upper limit R\ = 1, propagation is allowed 
only for sin 2 9 < |. As R\ decreases, the range of allowed propagation directions increases until it 
encompasses all angles at the lower limit R\ = ^. This polarization involves all three transitions 
and its growth rate is 



2(sin 2 9 + R l cos 2 9) + \ sin 2 29 J/J s 
'2 + (3 + sin 2 9)J/J S + sin 2 29 {J/J s ) 2 ' 



\ — K K] _ , ; "V J111 »-r-l.™° "J-r 4 ^ ^ "/ »* /c fi A 

^ — "7 — ^0 n /o • 2 /i\ t / t ■ 2 n/i c t / t \9 ' V d ~°J 



As with the polarization in a magnetic field, once the maser saturates for propagation at a given 
direction, the rate of further growth in the saturated domain is angle independent. 



5.2. Type Polarization 

Since this polarization involves all three transitions, R\ must adhere to the bound of eq. A12| . 



|1 — R\\ < Except for the absence of circular polarization, this solution is identical to that for 
overlapping Zeeman components (see equations 4-10| and 4-ll|) . The J-dependence of q can be 



neglected when J/J s > 2\1 — R%\, the resulting stationary polarization is 

2(3 - 2flx) 

q ~ l (2i? 1+ l)sin 2 0- (57) 

At .Ri = |, the upper limit allowed for this solution, it produces q = +1, same as the type 1 
solution that takes over at that value of R\. The requirement \q\ < 1 constrains the propagation 
directions to 

^^TP < 5 - 8 > 

the complementary angular range of the type 1 solution q = —1 (see eq. [5-5| and figure 3). At 
their common phase-space boundary, these two polarizations are the same. As with the analogous 
case of overlapping Zeeman components, this polarization is not stationary at low intensities. The 
requirement that \T n \ be smaller than 1 when r m = t s produces the constraint 

|(l- J R 1 )(sin 2 ^-|)|<-. (5-9) 

Figure 3 displays the bounds produced by this constraint for r s = 15 when they further restrict the 
phase space for physical solutions. Finally, the requirement that the source function be negligible 
when the J-dependent term dominates the behavior of K[ is 



J/J S Z \1-Ri\ > JS /J S , (5-10) 



the analogue of eq. 4-16 
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6. SUMMARY AND DISCUSSION 

The results presented here provide the general solutions for maser polarization in all cases that 
have been addressed in the literature thus far. Table 1 summarizes the properties of the solutions 
for any magnetic field strength, table 2 does the same for m-dependent pumping with any degree 
of anisotropy. These tabulations are further aided by figures 2 and 3 which display the phase space 
regions where each solution is applicable. Appendix B provides an analysis of the phase relations 
obeyed by the electric fields of the polarization solutions. 



6.1. Comparison With Previous Studies 

Previous detailed studies concentrated mostly on maser polarization in a magnetic field and 
were conducted in one of two limits, either xb ^ cxa or ig = 0. The results of §3 fully corroborate 
those of previous studies in the first limit (GKK, paper I). For any spin, this polarization is estab- 
lished at a given frequency once the Zeeman pattern is resolved at that frequency, the approach to 
the limit solution is at least exponential in xb- The polarization is established at the line centers of 
the Zeeman components when xb k; 1 and takes hold over the entire spectral region of interest when 
xb 2. A significant new result involves the absorption coefficients for 1 — * transitions (eq. 3-4 ). 
These show a considerable degree of anisotropy, introducing substantial differences between the 
behavior of the ir- and cr-components at most angles. 

Overlapping Zeeman components were previously studied only at the unphysical limit xb = 0. 
The xb — * limit of the solution derived here (§4.2) verifies the linear polarization found in 
the previous studies (GKK; papers I, II). However, while terms of higher order m xb have an 
insignificant effect on the magnitude of q, they produce some major new results. 



6.1.1. Circular Polarization 



The linear polarization is accompanied by circular polarization, proportional to xb (eq. ^j] 
12). For saturated masers, I follows the Gaussian frequency profile of the unsaturated absorption 
coefficients and the Stokes parameter V (= vl) is 

4T(0)*£ _ 
6 cos 



Here 1(0) is the intensity at line center and 



1 

x B 



,.2 



F(x;x B ) = — e~ {x+XB)2 - e"^-^) 2 ] ~ -4s e"** (6-2) 



is the frequency profile of the circular polarization. The last approximation is for xb < 1, corrections 
are only of order x\. Thus the profile is independent of xb to a high degree of accuracy in the entire 
region of interest, and is plotted in figure 4. The two peaks of the function occur at x = ±l/v2 
and the magnitude at the peaks is y/8/e. As the magnetic field strength varies, the profile of V is 
simply scaled by xb- The frequency separation of the two peaks, t/2A.vd, remains fixed and can 
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be used to determine the Doppler width, but not the field strength. Once the Doppler width is 
known, the Zeeman splitting can be determined from u pe ak) the ratio V/I at the peak of the Stokes 
parameter V, through 

X B = "^T «pcak COS 0. (6-3) 

As usual, because of the uncertainty about the direction of the field, its magnitude cannot be 
uniquely established. Since the spectra of the Stokes parameters obey V oc dl/dx, the magnetic 
field can also be determined using the method described by Troland & Heiles (1982) for data 
analysis when xb < 1- 

The circular polarization is responsible for an upper limit on the directions allowed for polarized 
maser radiation, and propagation perpendicular to the magnetic field is forbidden. The largest 
direction possible at Zeeman splitting xb and frequency x is 9 = cos~ 1 (4\^2xxb), the upper bound 
displayed in figure 2 for some frequencies x. Propagation along this boundary of allowed phase- 
space produces v = 0.94. Except for the immediate vicinity of this extreme edge, v is of order xb 
for most propagation directions. 



6.1.2. Onset of Linear Polarization 

In addition to introducing circular polarization, the higher order terms hold the key to some 
fundamental properties of maser polarization when xb < 1 that are not directly accessible to 
observations. At xb = 0, the J-dependence of the type solution for q disappears (see eq. 4-10| ; 
note that Ri = 1 in this limit), a singularity that caused considerable difficulty in previous studies 
because of the conceptual problem associated with polarization growth: Since the seed for maser 
radiation is spontaneous decays, when xb = the radiation starts unpolarized with intensity So. 
And since the polarization solution does not introduce any additional intensity scale in this limit, it 
was not clear at which intensity this solution would become applicable. Attempts to overcome this 
difficulty involved introduction of extraneous intensity scales through various arguments. While 
GKK claimed that the polarization sets in only after the maser saturates (J/J s > 1), paper II 
argued that J/J s p> 1/t s was the more appropriate condition (r s is the saturation optical depth). 

The general analysis presented here finally resolves this issue — the linear polarization solu- 
tion requires J/J s ^ xb > \J~Sq]~J s (eq. 4-16 ). The reason is quite simple. This polarization is 
established when the average populations of the m = ±1 and m = magnetic sub-states of the 
maser upper level become equal to each other. In the presence of Zeeman splitting and the absence 
of interaction with maser radiation, these populations are different. Interactions with properly 
polarized maser radiation equalize the sub-level populations once these interactions dominate the 
pumping process. This occurs when BJ/T (= J/J s ), the radiative rate relative to the loss rate 
of the pumping scheme, exceeds xb- By comparison, saturation is the process in which radiative 
interactions reduce the inversion itself, the difference between the average populations of the entire 
upper and lower levels, and requires J/J s > 1. The polarization solution requires amplification 
with an optical depth ^r s , halfway through the unsaturated domain. 

The intensity scale of saturated maser radiation is determined by J s ; the larger is J s the more 
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powerful the maser. For a given source function, eq. 4-16| shows that the polarization solution 
requires x B J s > Sq — stronger masers can enable the polarization solution at smaller Zeeman 
splitting. Conversely, when the maser pumping becomes weaker at a fixed Zeeman splitting, this 
condition cannot be met and the polarization disappears; as the maser begins to resemble a thermal 
source, it becomes unpolarized. 

6.1.3. Spurious Solutions 

While the xb — ► limit of the solution derived here reproduces the linear polarization of 
previous studies, there is one notable exception. None of the spurious solutions which required 
stability analysis for their removal in those studies is generated here. Indeed, as is evident from 
both tables 1 and 2 and figures 2 and 3, in any given region of phase space there is never more 
than one physical solution, if at all. In particular, the solution q = —1 whose stability properties 
for sin 2 9 < | caused some problems is not the xb — ► limit of any physical solution. Rather, it 
is only obtained as the xb — > limit of a spurious solution whose imaginary part vanishes in this 



unphysical limit (see §4.1, especially eq. 4-9). There are no physical solutions for polarized maser 
radiation at sin 2 9 < | when xb < 1, in agreement with the conclusions of paper II. 

It is worth noting that the results of a formalism that assumes xb = at the outset cannot be 
distinguished from those of the equally unphysical limit R\ = 1 of anisotropic pumping. Indeed, 
the latter case fully reproduces the GKK results. At sin 2 9 < | it has q = —1 and at sin 2 9 > | 
it has the same polarization as the xb = formal limit of the magnetic case. Obviously, though, 
anisotropic pumping has no direct relevance to the behavior in a magnetic field, as is evident when 
the two models are compared at physical, i.e., finite values of xb and |1 — R±\, however small. 
This further underscores the inherent pitfalls in a formulation that assumes a limit from the outset 
instead of properly approaching that limit within a general framework. 



6.2. Polarization Variation with Field Strength 

With the dependence on xb explicitly included, the evolution of the polarization with magnetic 
field strength can be described. The bounds placed on the Zeeman splittings and angles that allow 
different physical solutions are displayed in fig. 2 and table 1. At xb 10~ 3 , the allowed physical 
solution corresponds to the type polarization listed in eq. 4-12) , applicable for essentially all 



frequencies and sin 2 9 > |. As the field strength increases, the allowed angular region shrinks at 
a frequency-dependent rate that increases with distance from line center, thus the spectral region 
of allowed propagation shrinks inward and the line wings disappear first. This can lead to narrow 
maser features with linewidths less than the Doppler width in sources with xb ^ 0.1. Such features 
will display substantial circular polarization. Finally, when ^ 0.7 propagation is no longer 
possible for any direction even at line center, and the entire line is unpolarized. 



All along, the type 1 polarization listed in eq. 4-5 provides another allowed solution, albeit one 



that is unstable and that does not grow into the saturated regime so long as xb < 1. As the field 
strength increases further, these deficiencies are removed and this solution becomes the physical 
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polarization listed in equations 3-1 and 3-2, first at the centers of the three Zeeman components 
when xb > 1, spreading to the line wings when xb > 2. This x^-independent solution describes 
the physical polarization at all higher values of xb- 



6.3. Higher Spins 

Whenever explicit expressions for the absorption coefficients were required, those of J = 1 — ► 
transitions were utilized. Absorption coefficients for all other spins, which follow directly from 
general level population equations listed in paper I, can be avoided in most applications because 
various properties of the solutions are independent of the transition spins. In particular, type 
polarization is characterized by absorption coefficients that are equal to each other (eq. ^| 
H), and because of general properties of the vector-coupling coefficients, such solutions are spin 
independent (paper I). Thus the type solutions, in particular the polarization in a magnetic field 
with xb < 1 3 are valid for all spins (cf papers I, II). Further, the polarizations of fully resolved 
Zeeman patterns (equations 3-1 and |3-2f ) were derived without even specifying the absorption 
coefficients, thus they too are spin independent. Only the detailed expressions for the absorption 
coefficients of the ir- and er-components listed in equation 3-4 are based on explicit results for 
1 — ► transitions. These expressions are the starting point for analysis of the intensity angular 
distribution for the different Zeeman components, so this discussion strictly holds only when the 
1 — ► results are applicable. Included in this category are all transitions between equal spin states 
and m-independent pumping, because they share the same simple Zeeman pattern, so these results 
apply without modification to the OH main lines but not to its satellite lines. Still, the general trend 
of anisotropy is expected to hold in the latter case too (note in particular the angular dependence 
of the general expressions for the growth rates in equations |3-1| and |3-2| ) . 

It is important to emphasize the fundamental assumption that the only degeneracy of the 
transition levels involves their magnetic sub-states. When any of the maser levels includes additional 
degeneracy, so that magnetic sub-states of different levels overlap, the tight constraints responsible 
for the general solutions no longer apply and the polarization can be expected to disappear. This 
can be easily understood by considering the imaginary limit in which the hyperfine splitting of 
the OH molecule vanishes and its four ground-state maser lines blend into one. The polarization 
properties of masers with hyperfine degeneracy, notably H2O and methanol, cannot be studied 
without specifying the full details of their pumping schemes. The general solutions derived here 
apply to such masers only when they involve excitation of a single hyperfine transition. 



6.4. Comparison With Observations 

Since most results of the previous studies are reproduced here, the basic overall successes of 
maser polarization theory remain intact. OH is the only astronomical maser whose Zeeman pattern 
can be fully resolved with typical interstellar fields. It fulfills the condition xb > 1 for magnetic 
fields of order milligauss, frequently measured in HII/OH regions. Indeed, these masers are properly 
described by the theory for resolved Zeeman patterns. The most detailed polarization measurements 
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were done for W3(0H) by Garcia-Barreto et al. (1988) and Bloemhof, Reid & Moran (1992), and 
are adequately explained if Faraday rotation is responsible for removal of some linear polarization. 
A significant new result is the great disparity, which persists into the saturated domain, between 
Zeeman components at different angular regions, resulting in a preponderance of cr-components 
(§3). Garcia-Barreto et al. find that there are no features in W3(OH) that might be identified as 
7r-components, even accounting for possible Faraday rotation. The results of §3 adequately explain 
this observation if the magnetic field is aligned within less than 55° from the line of sight. 

For overlapping Zeeman components, the conclusions of papers I and II regarding linear po- 
larization are reproduced and the observational implications of these studies carry over. Maser 
polarization is established during the unsaturated exponential growth phase and is independent of 
spin. All pure-spin maser transitions of non-paramagnetic molecules should display the same polar- 
ization properties. This appears to be the case for SiO masers, which generally display high degrees 
of linear polarization. Mcintosh & Predmore (1991, 1993) have extended SiO linear polarization 
measurements up to J = 3 — ► 2 and their observations indicate that polarization properties indeed 
are J-independent — different lines display detailed similarities between fractional polarization, 
polarization position angle, and rotation of position angle with velocity. Faraday depolarization 
affects the low rotation states more because of their longer wavelengths and the polarization can be 
expected to decrease toward lower angular momenta, as observed. Although detailed calculations 
of Faraday rotation have yet to be performed, Mcintosh & Predmore (1993) find this to be the 
most plausible explanation of their data. 

Unlike SiO, which is a simple rotor, H2O maser transitions involve overlapping hyperfine 
components and the polarization can be expected to disappear in general. Indeed, H2O masers 
are generally only weakly polarized (e.g. Barvainis & Deguchi 1989). Exceptions do exist, though, 
and H2O masers sometime display high polarization, notably during outbursts such as the one in 
Orion (Garay, Moran & Haschick 1989; Abraham & Vilas Boas 1994). Such cases may involve the 
excitation of a single hyperfine component. Similar to H2O, methanol lines are in general split 
into many hyperfine components due to a number of mechanisms and can be expected to have 
polarization properties similar to H2O. Indeed, this is what Koo et al. (1988) found in a methanol 
12.2 GHz survey. This transition has a spin-rotation splitting of only ~ 2-3 kHz (Gaines, Casleton 
Sz Kukolich 1974), a degeneracy that could be responsible for the similarity to H2O polarization 
properties. 

A major new result is the circular polarization for xb < 1- The Stokes V spectrum has the 
traditional Zeeman anti-symmetric S-curve profile, plotted in figure 4. Such spectra are not always 
observed. However, this profile follows directly from the general symmetry in left-right rotations 
about the magnetic axis, independent of the specific theory developed. While this symmetry cannot 
be broken with the quiescent medium and uniform magnetic field assumed here, the model can be 
supplemented with symmetry breaking in the form of magnetic and/or velocity gradients, as done in 
various filter mechanisms (e.g. Deguchi & Watson 1986, and references therein). Such modifications 
can easily account for the removal of one component. 

Circular polarization, at the level of a few percent, was discovered in SiO observations of late- 
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type stars by Barvainis, Mcintosh & Predmore (1987). Following their discussion we take 9 ~ 45° 
and denote the percent maximal circular polarization by m c (= 100i> pea k), then eq. |6-3] shows that 
xb = 2xl0~ 3 m c . Observed values of m c range from 1.5 to 9, therefore xb varies from 3xl0 -3 to ~ 
10" 2 . The magnetic field can then be determined from eq. |1-1| , which leads to B = 2m c Av5 G for 
the 43 GHz line. Invoking a filter mechanism proposed by Deguchi & Watson (1986), Barvainis et 
al. devised an estimate for the magnetic field that has an identical expression, only the numerical 
coefficient is a factor of 8 larger .i Therefore, once their estimates for individual sources are scaled 
downward by an overall factor of 8 they hold as is, resulting in fields of 2-10 G. 

When the transition frequency varies, the Doppler width varies proportionately while the 
Zeeman splitting remains fixed. Therefore xb is inversely proportional to frequency (cf eq. 1-1) 
and the circular polarization decreases with the transition frequency when all other properties 
remain fixed. The circular polarization of SiO masers can be expected to decrease when the 
rotation quantum number increases. Indeed, observations of VY CMa by Mcintosh, Predmore & 
Patel (1994) show that the circular polarization of J = 2 — > 1 is smaller than that of J = 1 — ► 0, 
the opposite of the trend displayed by linear polarization. 

Circular polarization is observed also in OH maser emission from late-type stars. In the case 
of supergiants it can be substantial — 50% and higher is common (e.g. Cohen et al. 1987, and 
references therein). In the absence of theory for xb < 1, this was taken as a signature of xb <; 1 
and magnetic fields of at least a few milligauss. However, this polarization often appears as sharp 
reversals between adjacent narrow spectral components of the parameter V, similar to the profile 
displayed in figure 4, thus the two circular components are not separated and the data must be 
analyzed according to the theory for xb < 1- The separation between the peaks cannot be used as 
an indication of the magnetic field strength, neither is a determination of Doppler width from a fit 
to individual components an adequate procedure. Instead, the Doppler width should be determined 



from the separation of the two peaks, the field strength from eq. 6-3 



Since the OH circular polarization is typically an order of magnitude higher than for SiO, the 



parameter xb is generally an order of magnitude larger. Assuming for simplicity 9 ~ 45°, eq. 6-3 
shows that 50% circular polarization corresponds to xb — 0.1, consistent with overlapping Zeeman 
components. From figure 2, at such a large xb propagation at 45° is forbidden at x = 1 and is 
allowed only closer to line center. This illustrates the importance of constraints at such high values 
of xb and may help explain the narrowness of observed components in the profile of V. For the 
1665 MHz line, xb = 0.1 requires a field of only 0.4 milligauss for a Doppler width of 1 km s . 
For any polarization, the field can be obtained from B = 0.7u pea kAt>5 mG, and it is only ~ 0.1 mG 
for typical parameters. Thus the magnetic fields in supergiant maser regions are one to two orders 
of magnitudes lower than what was previously inferred in the literature (e.g., Reid et al. 1979, 
Claussen &: Fix 1982, Cohen et al. 1987). In contrast, maser emission from OH/IR stars displays 
somewhat lower circular polarization 15%), which was properly analyzed by Zell & Fix (1991) 
with the Troland & Heiles (1982) method. This analysis yielded estimates of field strengths < 



4 Note that the Barvainis et al. expression is written in terms of full width at half maximum. 
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0.1 G, and those remain valid. 

Circular polarization with i>p ea k ~ 10 -3 was measured for H2O masers in star forming regions 
by Fiebig & Gusten (1989). In general, the polarization solutions derived here do not apply to 
H2O when more than one hyperfine component is involved. Still, the profile of eq. |6-2| , displayed 
in figure 4, is a general one and can be expected to describe circular polarization whenever it is 
generated for xb <C 1. Fiebig & Gusten analyzed their data with the Troland & Heiles method, in 
essence relying only on the shape of this profile. Thus their estimates of magnetic fields of 15-50 
mG can be considered reliable. Analysis based on eq. |^3| produces similar results. 

The success of theory in explaining all the essential features of maser polarization observations 
is encouraging. Especially gratifying is the success in explaining the behavior of both circular 
and linear polarization for different maser transitions of the SiO molecule. These successes can 
be considered positive indication that all the relevant ingredients have been properly incorporated 
into the theory. A major component still missing is Faraday rotation. The study of this effect is 
left for future work. 

I thank R. Barvainis and K. Menten for useful information regarding SiO and methanol, 
respectively. This work was supported by NSF grant AST-9321847. 



A. APPENDIX: RADIATIVE TRANSFER AND ABSORPTION 

COEFFICIENTS 

General transfer equations for the Stokes parameters of line radiation were derived by Litvak 
(1975) and summarized in paper II. They are reproduced here for completeness. For any spin tran- 
sition, the radiative interactions with photons of polarization Am (= +1,0,-1) are characterized 
by an absorption coefficient K Am . Introduce 

k 1 = !(«+ + «-), (Al) 

the mean absorption coefficient for Am = ±1 transitions, and the three linear combinations 

Km = i[ K 1 (l + cos 2 (9) + K sin 2 6l] 
Ki = -K°)sin 2 6< 

k-c = — k~) cos9, (A2) 

where 6 is the angle of the radiation propagation direction from the quantization axis. Then the 
radiative transfer equations are 



dl 

dl 



K m I + KlQ + K C V 



= K mQ + KlI 



dU 
KT 
dV 
~dT 



KmY + K C I. (A3) 



-25- 



The basic absorption coefficients K Am are derived from the solution of the steady-state level 
population equations. As a concrete example, consider the standard model of a spin 1 — ► transition 
whose three upper levels are denoted (1, m) (so that Am = m). The loss rate V is assumed equal 
for all states, for simplicity, and the pump rates are P$ for the lower level and P% m (> Pq) for 
the (l,m) upper level. The pump rates are functions of the particle velocity distribution, assumed 
Gaussian. To simplify matters further, the geometry is assumed linear so velocities are directly 
related to frequency shifts from line center and the velocity Gaussian width translates to Ai/jj. 
Denote by J& m the angle- averaged intensity at the frequency of the (l,m) — * transition and by 
B the mean ^-coefficient. Then the steady-state populations, n\ m for an upper state and no for 
the lower state, are determined from 

P lm = Tn lm + Bf& m J Am (n lm - n ) m = Am = ±1, 

Po = Tn - fpJp( n ip ~ n o), (A4) 
p 

where 

/o = i(l-g)sin 2 #, /± = /i±ti; 

f x = ±(l + cos 2 6> + gsin 2 (9), w = \v cos B. (A5) 

The dimensionless factors f A m account for the radiative coupling of different polarizations (cf GKK, 
papers 1,11); note that /+ + /_ + / = 2/i + f = 1. 

The four rate equations provide the three equations for the population differences n\ m — no, 
which can be re- written as a set of equations for the three absorption coefficients K Am : 

(1 + fAmJA m /Js)K Am + E /Am'(4 m '/J S )« Am ' = ^ ■ (A6) 

Am' 

Here J s = T/sJl and (-Plm — -Po)/r are the appropriate unsaturated absorption coefficients, 

corresponding to the level populations in the limit J = 0. If we denote by a the sum on the 
left-hand-side, the absorption coefficients are 

Am _ 

K Am = ^ ° (A7) 

1 + J Am J A 

Inserting these results back into the summation defining a yields 

S CAm^™ , f AmJ Am /,„\ 

° = -, , ^ , where e Am = T , f 7 — • (A8) 

1 + e Am Js + J Am J Am 

This completes the definition of K Am in terms of the unsaturated absorption coefficients K Am and 
the angle- aver aged intensity. The results can also be expressed as rational functions 

K Am = N Am/ D ( A9 ) 



5 In the scalar maser model the definition is J s = T/(2B). 



-26- 



where 

D = 1 + 2 (jo + j+ + J- ) + 3 (jo 3+ + 3+3- + 3-3o ) + *joj+j- 
N+ = (1 + 2j + 2j_ + 3j j-)n+ - j_(l + jo)«o - io(l + J- 

^~ = -i+(l + Jo)4 + (1 + 2j + 2j+ + 3joi+)^o " io(l + 
N° = +(l + 2j+ + 2j_ + 3j+j_)«8 (AlO) 

and where j^m = f& m J&m/Js- 

The absorption coefficients of a polarized maser resemble those of three scalar masers whose 
unsaturated absorption coefficients are reduced by <r, an intensity-dependent common factor re- 
flecting coupling through the shared lower level. When f Am J Am <C Js (the maser is unsaturated at 
the frequencies of all three transitions), e^m — and a ~ so that K Am ~ /^ m , as should be. In 
the opposite limit, in which all three transitions are strongly saturated, e/\ m ~ 1 and a ~ | J] k^" 1 
so that 

Am - TJ^—^o" 1 ~ I E -o Am ')- (AH) 

"ijAmJAm Am'^Am 



K 



In order for all three transitions to remain inverted in the limit of strong saturation, the unsaturated 
absorption coefficient of each transition must exceed the sum of the other two by at least \ . From 
this it follows that the bounds 

Kg < + < 3«g (A12) 

must be obeyed whenever all three transitions couple at the same frequency, placing limitations on 
potential pumping schemes. 



B. APPENDIX: ELECTRIC FIELDS AND PHASES 

Maser photons are generated in stimulated emissions. An induced photon has both the fre- 
quency and wave vector of the parent photon, but not necessarily its polarization. Indeed, the 
polarization of the induced photon is determined by the change in magnetic quantum number of 
the interacting particle, not the polarization of the parent photon. Consider for example a photon 
produced in a Am = transition, and thus linearly polarized. Because such a photon can also 
be described as a coherent mixture of two circularly polarized photons, it can induce transitions 
with any value of Am, thus any polarization. For instance, in the case of a spin 1 — >• induced 
transition, when the interacting particle is in the m = state, the induced photon, too, is linearly 
polarized. But when that particle is in one of the m = ±1 states, the induced photon is circularly 
polarized, although the interaction amplitude is reduced. 

The polarization solutions correspond to waves with electric fields whose phases are properly 
tuned to the level populations determined by the pump rates so that induced radiation maintains 
the polarization of the original one. Methods to determine such electric fields were introduced in 
paper I and are easily generalized to the solutions derived here. Two coordinate frames are defined 
by the problem. One, the fc-frame, has its z axis aligned with the direction of wave propagation. 



-27- 



Components in this frame are denoted by subscripts. The other, the .B-frame, has its z axis aligned 
with the quantization axis and components in this frame are denoted by superscripts. The common 
y axis is perpendicular to the plane that contains the z axes of the two frames. Since the electric 
field of the wave is a vector, its components in either frame can be obtained from those in the 
other by a simple transformation. In the case of fully resolved Zeeman patterns, only one Am 
transition is coupled and only one of the S-frame components of the electric field does not vanish 
at a given frequency. There is no phase relation to be determined in this case and the solution is 
fully prescribed in paper I. 

When all three components overlap, as is the case for small is or anisotropic pumping, denote 
by E Am the proper circular component of the electric field in the i?-frame. The corresponding 
intensity I Am oc \E Am \ 2 obeys I Am = f/± m I. Therefore, if we introduce the complex unit vector 
e = E/|E|, its circular components in the 5-frame are 

=° = , FTn p+ = , /77p^+ p~ = . I f p i( t> 



e u = V/o, e + = ^/+e^+, e~ = yj f.e*- , (Bl) 

where 4>± are phase angles. The choice of a real phase for e° is always possible because the overall 
phase of e is irrelevant. With these definitions, the rectangular /c-frame components of e are 



e x = _ (y/ +e »*+ + ^/_ e *»- ) cos + V/o sinfl 

e z = ^(yfUe** + JJlj*-) sine + ^/ToCOs9. (B2) 

Propagating waves are transverse, i.e., e z = 0, and the normalized Stokes parameters are determined 
from the two non-vanishing components via q = e^ — e^, u = 2Re(ea;e*), v = 2Im(e x e*). Since the 
transverse condition e z = is a complex relation, it provides two equations. From the imaginary 
part, 

- ^JJl$m4> + = ^JJlsvcKp- = V7i sin0 (B3) 

so the field structure is determined by a single phase, 4>, whose magnitude is obtained from the real 
part of the transverse condition: 

2 _ tf + w 2 tan 4 1 - q 2 - v 2 

C ° S ^ 2/o/itan 2 (l-q)(l + cos 2 e + qsm 2 ey [ ' 

This completes the definition of the vector e for a polarization mode characterized by a specific 
set of q, u and v or, equivalently, /+, /_ and /q. The results are essentially kinematic relations 
obeyed by all polarization modes. Stationary modes are distinguished as those that maintain a 
constant unit vector e, requiring special phases. Type 1 polarizations have q 2 + v 2 = 1, therefore 
they have = 0, except for q = +1 which has 4> = i^ 71 " ( so this is never the v — > limit 
of another type 1 polarization). In the case of type solutions, the level population equations 
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determine the stationary, J-independent part of the structure functions as /o = (3 — 2Ri ) / (2R\ + 1) 
h = (2Rx - 1)/(2R 1 + 1). Therefore 



2Ri 



cos 



+ 



8R 2 C tan 2 9 



2(2i?i - l)tan 2 (9 (3 - 2 J Ri)(2i?i - 1) 



(B5) 



for all type polarizations, where R\ and R c are the ratios appropriate for the particular case. 
Various bounds and properties of the type solutions for small xb and for anisotropic pumping find 
simple explanations in terms of the variation of the phase 4> with pump rates and the requirement 
that it be real (i.e., < cos 2 </><!) for physical solutions. 



Table 1. Polarization Solutions — Magnetic Field 



Validity Domain 



Type 



y/S^/T s <x B < 1/V2; 
J I J s > x B ; 
sin 2 > ±[l + fx|(l + 2x 2 )]; 
cos 9 > \\[2xx B a 



1 + \x\{\ - 2x 2 



3 sin 2 (9 



16xx B 
3 cos 9 



1 + PJJs 



XB > 1 + X] 7T 



xb > 1 + x; <r~ 



sin 



1 + cos 2 





2 cos 9 
1 + cos 2 ' 



K sin 



2 a c 



\k{1 + cos 2 6») 



a These are the primary constraints on the 9—xb phase space for this solution (eqs. [4-3j , |4-13j and 
4-16). Additional constraints are provided by eq. |4-15| . 

b These results hold for all spins, with Kq the absorption coefficient in the absence of both Zeeman 
splitting and maser radiative interaction. 

c These results hold for all spins, with k the absorption coefficient of the transition. Detailed 
expressions displaying the explicit dependence of n on J/J s are listed for 1—^0 transitions in 
eq. 



3-4 



Note. — Polarization solutions for masers in magnetic fields characterized by Zeeman splitting xb 
(eq. |l-l|) . The different types of polarization solutions are discussed in §2.1.1. 
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Table 2. Polarization Solutions — Anisotropic Pumping 



Validity Domain Type 



\ < Ri < 1; 2(sin 2 9 + Ri cos 2 g) + | sin 2 29J/J S 

sm 2 6< Ift^+l ~ K ° 2 + (3 + sin 2 6)J/J S + sin 2 26{J/J s f 

J/J s > 2)1-^1; o 1- 2(3~2fli) «8(2fli + l)/3 

•2^3-2^1 (2 J R 1 + l)sin 2 1 + § J/J s 

sm y > 2RT+T 



Note. - Polarization solutions for masers in the case of pumping with a degree of anisotropy 
characterized by R\ (eq. |5-1[) . Kq is the unsaturated absorption coefficient of the Am = transition, 
k,q that of |Am| = 1. 
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FIGURE CAPTIONS 



Figure 1: The polarization vector II and the vector n p that controls its radiative transfer in the 
space defined by the normalized Stokes parameters q, u and v. 

Figure 2: Phase space domains of physical solutions in a magnetic field with Zeeman splitting xb 



(eq. |l-lD . At each indicated frequency shift x = (v — Uq)/Avd : polarized maser propagation is 
allowed only inside the 6-xb region enclosed by the corresponding boundary. The domain for each 
x contains all the domains for larger values of x. The primary bounds when xb < 1 are from 



eq. 4-13 and are denoted with full lines. At some frequencies they are supplemented by the bounds 
of eq. [4-3| (long-dashed line) and eq. |4-15| (short-dashed lines), assuming r s = 15. An absolute lower 
limit on xb for all frequencies, dependent on the pumping scheme and not plotted, is provided by 



eq. 4-16. When xb > 1, propagation is allowed for xb > 1 + x and all directions. The thin dotted 



line in this region marks the separation between directions where ir- and cr-components dominate. 
Figure 3: Phase space domains of physical solutions for m-dependent pumping whose degree of 



anisotropy is characterized by R\ (eq. 5jl). The solutions applicable in the different regions are 
listed in table 2. The short-dashed lines that further restrict the corners of the domains of type 
polarization are from eq. |5-9| , assuming t s = 15. The dotted lines on the two sides of R\ = 1 are 
a schematic representation of the bound of eq. |5-10 , The actual location of these lines depends on 
the pumping scheme. 

Figure 4: Profile of the circular polarization when ig <C 1 (eq. |6-2j). 



